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Recent developments is the field of robust multivariable control have merged the theories of and Hi control* 

This mixed Hz/Hoo compensator formulation allows design for nominal performance by ff 2 norm minimization 
while guaranteeing robust stability to unstr u ct ur ed uncertainties by constraining the Hoc norm. A key difficulty 
associated with mixed Hz/Hoc compensation is compensator synthesis. A homotopy algorithm is presented for 
synthesis of fixed order mixed Hz/Hoc compensators. Numerical results are presented for a four di«k flexible 
structure to evaluate the efficiency of the algorithm. 


Introduction 

M ODERN control theory has revolutionized control system 
design, with H 2 and methods gaining widespread recog- 
nition and application in controller synthesis for single-input/single- 
output (SISO) and multi-input/multi-output (MIMO) problems. 
Early work in multivariable control synthesis utilized a quadratic 
cost functional to minimize the 2 norm of a system response to white 
noise inputs. Although the H 2 procedure is well suited to many sys- 
tems that specify performance in terms of rms quantities such as 
minimizing line-of-sight errors or control energy, it is well known 
that stability and performance cannot be guaranteed in the presence 
of model uncertainties. Robustness is addressed in control the- 
ory, which guarantees stability and performance (when defined by 
an co-norm measure) in the presence of unstructured uncertainty 
models, albeit often resulting in overly conservative designs. A sig- 
nificant disadvantage of these modem control techniques is that the 
resulting compensator is the same order as the generalized plant, 
which is often larger than the original plant due to the inclusion of 
frequency dependent weights to achieve the desired performance 
and robustness characteristics. 

The consequential large controller order can be indirectly allevi- 
ated by reducing the order of the controller or alternatively by re- 
ducing the order of the design plant In either case, indirect methods 
are suboptimal in performance and do not guarantee closed-loop 
stability. However, direct methods may be employed that impose 
constraints on controller order or architecture in the optimization 
procedure and, hence, provide stability and performance guarantees. 
In an optimization-based synthesis procedure, necessary conditions 
are formulated for the constrained closed-loop system that ensure 
internal stability. The optimal projection approach 1 is an Hi proce- 
dure whereby order constraints are imposed on the controller and 
the necessary conditions for minimizing a quadratic cost functional 
with respect to the fixed-order controller are derived. The resulting 
necessary conditions consist of two modified Riccati equations and 
two modified Lyapunov equations coupled by an oblique projection 
matrix. However, solution of the necessary conditions for realistic 
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large-order systems is a difficult task. Homotopy methods have been 
employed to solve the optimal projection equations. 2 

As a means of providing robustness and performance, the mixed 
Hi/ Hoc methodology has been developed. Much work has been 
done with variations of the mixed H 2 f Hoc problem. (For a sum- 
mary, see Ref. 3.) The earliest refereed work was done by Bernstein 
and Haddad, 4 who extended the optimal projection approach to lin- 
ear quadratic Gaussian (LQG) control with an norm constraint. 
Their formulation minimized an overbound on the H 2 norm from a 
disturbance input to one output while satisfying an Hoc norm over- 
bound from the same disturbance input to a second output. The 
difficulty with this approach is the potential conservatism resulting 
from minimizing an overbound on the H 2 norm. The first attempt 
at solving the general mixed H 2 / problem was by Rotea and 
Khargonekar, 5 who allowed independent inputs and outputs for the 
two transfer functions and minimized the actual H 2 norm based 
on full state feedback. Ridgely and Walker 3 extended the formula- 
tion to output feedback including the fixed-order case with either 
regular or singular constraints, and Canfield et at. 6 provided a 
numerical solution. Another approach to the general mixed H 2 j Hoc 
problem was developed by Sweriduk and Calise, 7 who used a dif- 
ferential games formulation to obtain fixed-order controllers. A nu- 
merical solution of this formulation using homotopy was developed 
by Whorton et al., 8 an extension of which is the subject of this paper. 
Davis et al. 9 also recently developed an algorithm for mixed H 2 / Hoc 
design. 

Ridgely and Walker 3 showed the mixed H 2 /H^ to be a strictly 
convex optimization problem with a unique solution when the con- 
troller is of the order equal to or larger than the underlying H 2 
problem. Tlie solution is shown to lie on the boundary of the in- 
finity norm constraint when active (for y < y, where y is the 
norm when the optimal H 2 controller is used) and is just the H 2 
controller when y > y. For controllers with order less than the 
underlying H 2 problem, the possibility of local minima in the un- 
constrained H 2 problem (optimal projection) exists and the solution 
of the mixed H 2 /H^ problem may or may not lie on the boundary 
of the Hoc norm constraint As will be shown in the next section, 
Sweriduk and Calise 7 begin with the fixed-order Hoc cost functional 
and append the fixed-order H 2 cost functional. Reference 3 takes the 
converse approach by appending the norm constraint to the H 2 
cost functional. As a consequence, the formulation in Ref. 3 can han- 
dle both regular and singular //„ constraints and can be specialized 
to the H 2 problem. Conversely, the Ref. 7 formulation can handle 
both regular and singular H 2 constraints and can be specialized to 
the Hoc problem. However, the ability to handle singular constraints 
on either cost exists as long as one of the cost constraints is regular. 
Whereas the Ref. 7 formulation assumes that the plant dynamics 
( A matrix) of the two transfer functions are the same. Ref. 3 al- 
lows different dynamics in the two separate underlying //«» and H 2 
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problems. Using the canonical compensator with a static gain output 
formulation presented in following section, the Ref. 7 formulation 
requires the simultaneous solution of five coupled nonlinear matrix 
equations. Reference 3 presents the necessary conditions in the form 
of seven coupled, nonlinear matrix equations. 

The objective is to build on the results of Ref. 7 by presenting a 
homotopy algorithm that solves the mixed H 2 /H z c fixed-order com- 
pensator synthesis problem. Fixed-order H 2 and H ^ controllers are 
obtained as special cases of the algorithm. The paper is organized 
as follows. First, a formulation of the problem with the compen- 
sator in controller canonical form is presented and the necessary 
conditions for the fixed-order H 2 controller are developed. These 
results are then extended to the fixed-order H 0 0 and mixed U 2 j H x 
controller design using the differential game results of Ref. 7. (The 
approach and algorithms of this paper also follow Refs. 10 and 1 1 
for H 2 controller synthesis.) Second, homotopy methods are intro- 
duced and a homotopy algorithm is developed to synthesize mixed 
Hz/Hoc compensators. Finally, numerical results are presented for 
evaluation of these homotopy algorithms followed by a discussion 
and concluding remarks. 

Problem Formulation 

The generalized plant of a standard control problem is given by 


where the disturbances are confined to the set of signals with 
bounded power and fixed spectra. If the disturbance is modeled 
as white noise, the objective is 

min | J(A C , B c , C c ) = lim £{z(f) r z(f)} } (14) 

It can be shown that the cost can be expressed as 

J{A C , B c , C c ) = tr (QBB t ) = tr{PC r C} (15) 

where 

AP + PA T + BB t = 0 (16) 

A t Q + QA + C t C = 0 (17) 

P is the controllability jrammian of (A, B), and Q is the observabil- 
ity grammian of (C, A). An equivalent cost functional also arises 
for the case of impulsive inputs. 

To obtain the H 2 optimal compensator, the Lagrangian is defined 
as 

C(Q, L, A c . B c , C c ) = tr {QBB T + (A 7 Q + QA + C T C)L } 

(18) 


x = Ax + -i- B 2 u 

(i) 

z = C { x + D n u 

(2) 

y — C 2 x + D 2 iW + £> 22 H 

(3) 

where x e R n is the state vector, w € R nw is the disturbance vector, 
u € R nu is the control vector, z € R nz is the performance vector, 
andy € R ny is the measurement vector. The following is assumed. 

1) (A, 2?i, Ci) is stabilizable and detectable. 

2) (A, B p , C p ) is stabilizable and detectable. 

3) D\ 2 has full column rank. 

4) D 2 \ has full row rank. 

A general compensator for this system is 



x c = A c x c + B<y 

(4) 

U = CcX c 

(5) 


where x c e R nc is the state vector of the controller the dimension 
of which can be specified. Closing the loop using negative feedback 
yields the closed-loop system dynamics 


where L is a symmetric matrix of multipliers. Matrix gradients are 
taken to determine the first-order necessary conditions, 



dC 

dB c 


9 £ SC_ 

dL ~ ’ a A e 


0 , 



(19) 


Hence, computing an H 2 optimal controller of fixed order nc < n 
for the general controller structure given in Eqs. (4) and (5) requires 
the simultaneous solution of five coupled equations. This is not only 
computationally expensive but is also further complicated by the fact 
that the problem is overparametrized with such a compensator. 

A controller form architecture 12 is imposed on the compensator 
dynamics. This minimal realization avoids the problem of over- 
parametrization and is a canonical form under mild conditions. 13 
The internal structure of the compensator is prespecified by assign- 
ing a set of feedback invariant indices v, . In controller canonical 
form the compensator is defined as 


x c = P°x c + N V - N°y (20) 


x = Ax + Bw 
Z = Cx 


( 6 ) 

(7) 


u c = -Px c 

(21) 

u = —Hx c 

(22) 


where 




C = [C, -Z>| 2 CJ (11) 

The set of all internally stabilizing compensators is defined as 

S c = {(A c , B CJ C c ) : A is asymptotically stable) (12) 


where x c € R nc andw r € R ny . P and H are ffee-parameter matrices, 
and P° and N° are fixed matrices of zeros and ones determined by 
the choice of controllability indices. Similarly, a compensator in ob- 
server canonical form can be constructed. Only the controller canon- 
ical form is employed, which imposes the lower bound nc > ny on 
the order of the compensator. 

Let 

*=[“,] (2s> 

The augmented system may be expressed as 

(24) 

= Ax + B\W + B 2 u 


For an H 2 problem, the objective is to minimize the H 2 norm on 
the closed-loop transfer function from disturbance inputs to perfor- 
mance outputs 


T zw = C{sI - AY'B (13) 


Z = [C, 0]f + [Dn 0]« 

= C\i+D n u 


( 25 ) 
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y = 10 I]X = C 2 X (26) 

= -Gy (27) 

Equations (24-27) define a static gain output feedback problem 
where the compensator is represented by a minimal number of free 
parameters in the design matrix G. The closed-loop system is given 
by 



x — (A — B 2 GC 2 )x *f* B\W 
= Ax + Bw 


(28) 


z = (Ci - D l2 GC 2 )x = Cx (29) 

Minimizing the H 2 norm of T zw - C(sl - A)~ l B utilizes the same 
Lagrangian as given in Eq. (18), but now £ is only a function of 
three parameter matrices, i.e., C(Q , L, G). Thus, only three first- 
order necessary conditions result: 

3£ - -r --t 

— = AL + LA t + BB t = 0 (30) 


3£ 

— - = A r Q + QA + C T C = 0 
aL 


(31) 


3C 


— = 2 (Dj 2 D u GCi - Dl 2 C, - B[ Q)LC[ = 0 (32) 


t 


W 


u 



A 

c 

B 

e 



C 

c 

0 



Fig. 1 Mixed problem. 
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where 


r b p ' 


B, ~[ 

(42) 

II 

© 

(43) 

D lp = lD, p 0] 

(44) 


Controller canonical forms can also be used to solve the //«, 
problem. The objective is now to minimize the oo norm of the 
transfer function from disturbance inputs w to performance outputs 
z given in Eq. (13). In this case the necessary conditions for an 
suboptimal fixed-order compensator gain G are 7 

— = (A+,-WS.)l 

+ L(a + y~ 2 BB t Q 00 ) T + BB t =0 (33) 

BjC 

— =A T Q x + Q x A + C T C + y- 2 Q oo BB T Q x , = 0 (34) 
3£ 

— = 2(OfjD 12 GC 2 - D[ 2 C, - Bl Q x )LC t 2 = 0 (35) 

where 

£(G<*>, L, G) = tr [Q x BB t + ( A t Q x + Q X A 

+ C T C + y~ 2 Q x BB T Q x )L} (36) 

As in the H 2 problem, three coupled equations have to be solved 
to obtain a fixed-order compensator that satisfies the constraint 
II^zujIIoo < y • 

Fixed-order //«> design has also been extended to fixed-order /x 
synthesis. 1415 Because controller design is a subproblem when 
designing for robust performance with structured uncertainty, the 
fixed-order technique just introduced has the potential to constrain 
the order of the controller that is normally subject to significant 
increases in the /x- synthesis procedure. 

The mixed // 2 /ffoo problem can be approached in a similar 
fashion. 7 In this case, the generalized plant has additional inputs 
and outputs w p and z p , respectively, which define the H 2 perfor- 
mance criterion. The plant dynamics are illustrated in Fig. 1. The 
inputs w and outputs z are used to define the //«*> performance crite- 
rion. Using the controller canonical form for the compensator, the 
augmented system for the mixed problem is 


Ax + Biw + B p w p + B 2 u 

(37) 

z p — Cpi + Di p u 

(38) 

z — C\x + D\ 2 u 

(39) 

y = c 2 x 

(40) 

u — —Gy 

(41) 


The other expressions are the same as in Eqs. (24-27). Consequently, 
the closed-loop system is given by 

x = (A - B 2 GC 2 )x + B p w p + B\W 

_ . _ (45) 

— Ax- 1- BpW p + Bw 

z p = (C p - D lp GC 2 )x = Cpi (46) 

Z = (Ci - D l2 GC 2 )x = Cx (47) 


To formulate the performance index of the mixed problem, the 
Lagrangian for the H 2 problem in Eq. (18) is adjoined to the 
Lagrangian for the problem in Eq. (36) by a scalar weight X: 


JO = ti{Q x BB T + ( A t Q x + Q X A 


+ C T C + y- 2 QooBB T Q x )L + kXC T p C p 
+ (AX + xa t + B p Bj)L p } 

(48) 

The weight X on the H 2 norm allows a tradeoff between performance 
(H 2 norm) and robustness (H^ norm). The first-order necessary 
conditions are 

< f- = (A + y- 2 BB T Q x )L 

0 Stfoo 


+ L(A + y- 2 BB T Q x ) T + BB T = 0 

(49) 

3£ 

— = A r Q x + Q X A + C T C + y~ 2 Q x BB T Q x 

= 0 (50) 

dC - T - 

— = A t L p + L p A + \C T p C p = 0 

(51) 

3£ 

— = AX + XA T + B p B p = 0 

(52) 


— = l[bJ 2 D n GC 2 LCl - D\ 2 C x LC t 2 - B t 2 Q x LC t 2 

+ kb] p b ip GC 2 XC T 2 - kb T lp C p XCl - B T 2 LpXC T 2 ] = 0 

(53) 

Promising results have been obtained for the and the mixed 
problem where a conjugate gradient method was used in the 
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computation. 7 A disadvantage of this method is that convergence 
slows down near the optimum. Also, an initial starting guess for the 
compensator gain G has to be provided that stabilizes the closed- 
loop system. In this paper, a homotopy method is used to continu- 
ously deform the solution of a simple problem formulation to the 
solution of the desired problem formulation. 

Homotopy Methods 

Homotopy methods offer an attractive alternative to more stan- 
dard approaches of optimal controller synthesis such as sequential 
and conjugate gradient methods. The basic philosophy of homotopy 
methods is to deform a problem that is relatively easily solved into 
the problem for which a solution is desired. 

Homotopy (or continuation) methods, arising from algebraic and 
differential topology, embed a given problem in a parametrized fam- 
ily of problems. More specifically, consider sets U and Y e 31” and 
a mapping F: U — ► T, where solutions of the problem 

F(u) = 0 (54) 

are desired with u € U and F(u)eY. The homotopy function is 
defined by the mapping H:Ux[ 0, 1] -> 9^" such that 

(Ki, 1) = F(u) (55) 

and there exists a known (or easily calculated) solution u 0 , such that 

H(u 0 , 0) = 0 (56) 

The homotopy function is a continuously differentiable function 
given by 

H[u(ct), a] = 0, Va € [0, l] (57) 

Thus, the homotopy begins with a simple problem with a known 
solution, Eq. (56), which is deformed by continuously varying the 
parameter until die solution of the original problem, Eq. (54), is 
obtained. 16 The power of homotopy methods is that minimization 
is not strongly dependent on the starting solution but depends on 
local, small variations in the solution. Theoretically, these methods 
are globally convergent for a wide range of complex optimization 
problems, but in actuality, finite wordlength computation often in- 
troduces numerical ill-conditioning, resulting in difficulties with 
convergence. In light of these numerical limitations, a judicious 
choice of the initial problem and the associated initial stabilizing 
compensator is necessary for convergence and efficient computa- 
tion. However, the ability to select an initial problem with a simple 
solution renders homotopy methods more widely applicable than 
sequential or gradient-based methods, which have a stringent re- 
quirement for an initial stabilizing solution. 

Both discrete and continuous methods are used to solve the homo- 
topy. Discrete methods simply partition the interval [0, 1] to obtain 
a finite chain of problems 

H(u , a„) = 0, 0 = of 0 < <*i < • * • < (*n ~ 1 (58) 

Starting with a known solution at a n , the solution for H(u , a„ + 1 ) is 
computed by a local iteration scheme. Continuous methods involve 
integration of Davidenko’s differential equation, which is obtained 
by differentiating Eq. (57) with respect to a, yielding 



Given «(0) = this initial value problem may be numerically in- 
tegrated to obtain the solution at a = 1 if the solution exists and is 
uniquely defined. 

Research remains to be done in the application of homotopy algo- 
rithms. Efficient application of homotopy methods depends on the 
initial problem, the final problem, and the deformation undertaken. 
Given a good initial solution, the key to convergence is the ability to 
accurately track the solution curve, which is determined by the de- 
formation undertaken. The ability to predict the solution along the 


homotopy path via Davidenko’s differential equation makes con- 
tinuous homotopy methods superior to discrete methods. Efficient 
computation of die Hessian is the primary issue for practical im- 
plementation of continuous homotopy algorithms. In the following 
sections, a continuous homotopy algorithm is presented for fixed- 
order mixed H 2 /H oo compensator design. The set of homotopy al- 
gorithms for synthesizing fixed-order // 2 , H x , and mixed Hi/H^ 
compensators has been organized into a MATLAB™ toolbox called 
fixed-order compensation of uncertain systems (FOCUS). 

Homotopy Algorithm 

This section describes the algorithm used for implementing the 
continuous homotopy. In essence, a mixed discrete and continu- 
ous approach is employed where Davidenko’s differential equation, 
Eq. (59), is integrated along the homotopy path, and at discrete points 
along the trajectory, a local optimization is used to remove integra- 
tion error. The optimization scheme is a hybrid Newton/conjugate 
direction search method developed for this application. 17 Local op- 
timization at discrete points along the homotopy trajectory allows 
a crude integration procedure with large step sizes to be employed 
for efficiently tracking the solution curve. This approach follows 
closely that of Refs. 10 and 1 1 and is employed in the following 
algorithm. 

1) Find initial solution (a = 0). 

2) Advance a 

+ &<xo,k 

3) Predict 9 

9(oc\ k ) — 0( a 0 ) + Ao£ O> *0'(ao) 

where 



4) Check prediction error 

e k (6 t <x)=\\Je[0(ot Uk m < € 

a) if error less than tolerance, continue; b) if not, 0.5 Aao.* -► 
Aa 0 , *+ 1 ; and c) increment k and repeat steps 2-4. 

5) Correct with hybrid Newton/conjugate direction method to 
compute local minimum. 

6) If <x — 1, stop. Otherwise, go to step 2. 

Various approaches may be taken when selecting the deformation, 
but the general procedure applied in this effort is outlined as follows. 

1) Synthesize a low-authority H 2 (full-order) compensator. 

2) Reduce the compensator to desired order and transform to 
canonical form. 12 

3) Set y large enough so that the H 2 and compensators are 
approximately equivalent. 

4) Use homotopy to deform the initial low -authority, reduced- 
order Hi compensator into a full-authority reduced-order H 2 com- 
pensator (Hi homotopy). 

5) Deform the full-authority H 2 compensator into a nearly optimal 
Hoc compensator with y approaching its infimum (H^ homotopy). 

6) At discrete values of A., fix y and deform the compensator 
into the mixed H 2 /Hoo compensator by varying X (mixed H 2 /H x 
homotopy). 

This procedure was chosen because it has been observed numer- 
ically that order reduction techniques tend to work best for low- 
authority LQG controllers. 10 * 18 A canonical form is imposed on the 
compensator structure to minimize the number of free parameters, 
which in some cases can also lead to numerical ill conditioning. A 
balancing transformation, which does not affect the controller char- 
acteristics, relaxes the strict structure in the P° and N° matrices in 
Eq. (20) and improves the conditioning of the problem. 

The procedure outlined separates the compensator synthesis into 
distinct phases. The initial reduced-order full -authority compensator 
is synthesized using the H 2 homotopy, which is then deformed into 
the reduced-order H 0 0 compensator. During the H^ phase, the scalar 
H 2 norm weight X is fixed (as are the plant matrices) and only the 
Hoc norm overbound y is varied. At discrete values, y is fixed and X 
is varied to perform the H 2 norm minimization. Thus, the procedure 
alternates between the H ^ and H 2 norm minimization. 
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During the homotopy, both the predicted and corrected gains are 
checked to ensure closed-loop stability. After each correction step, 
the cost gradient is checked to verify descent During the ho- 
motopy, the solvability of the Riccati equation using predicted or 
corrected gains must also be checked. If any of these conditions 
are violated during correction, the correction step size is scaled and 
the condition is checked again. If scaling the correction step size 
is ineffective, the prediction step size is decreased and the predic- 
tion phase is repeated. This process continues until the homotopy 
is completed or until the prediction step size is decreased below a 
prespecified tolerance. 

The following sections detail the derivations employed in the 
homotopy algorithm for mixed H 2 / //«, design. A complete devel- 
opment of the Hi, Hoc, and mixed Hj/H^ algorithms is given in 
Ref 8. 


The derivatives with respect to 9 are denoted 


<*>“>= f 


( 68 ) 


For 6, = g ik , 


G 0) = ^ = £,* (69) 

which is a matrix of zeros except for a one in the ik element. 

To obtain expressions for L y) , Q$, L p \ and X (J) , differentiate 
Eqs. (49-52) with respect to 6j to obtain 

0 = (A + y- 2 BB T Q txl )L u) + L U) (A + y- 2 BB T Q x ) T 


Mixed Hz/Hoo Development 

The homotopy function as well as the gradient and Hessian matri- 
ces are determined from the first-order necessary conditions for an 
optimal mixed Hi/H^ compensator given by Eqs. (49-53). Define 
9 to be a vector comprising the free compensator parameters 

9 = vec(G) (60) 

where G is the output feedback gain matrix defined in Eq. (27). The 
gradient of the cost is 


+ [( A °> + y - 2 bb t q<£)l + l(a + y^bbTq'sY] 

T (70) 

0 * {a + y- 2 BB T Q x ) Q<£ + Q<£(A + y~ 2 BB T Q oa ) 

+ [A^Qoo + QocA ij) + (C r C)0>] (71) 

d-A T L^ + L^A + [A^ T L p + L P A U) + X(CjC p ) 0> ] (72) 

0= A X u) + X U) A t + [A 0) X + XA U)T + (B,Bj) 0) ] (73) 





= 0 


where 3£/3G is given by Eq. (53). 

The homotopy function is defined as 


H(9,a) = 


d£(9 , a) 

3 9 



(61) 


(62) 


Derivatives of the closed-loop matrices are obtained from Eqs. 
(45-47) and are given by 

A (j) = -B 2 C {j) C 2 (74) 

(BB t ) U) = (B p B T p ) U) = 0 (75) 

(C r C)°> = (ClG u)T Dj 2 D l2 GC 2 ) + (C^G u)T Dj 2 D n GC 2 ) T 


Note that L is now a function of the homotopy parameter or since 
the system matrices are now functions of a. The gradient of the 
homotopy function is 

V[ff r (0,cr)] = [V,// r V 0 // r ] (63) 

Computation of Hessian 

The derivative of the N x 1 vector valued homotopy function, 
H t (9) = [hi(9),h 2 (9 ) f . . . , h N {9)], with respect to the N param- 
eter vector 9 is the N x N Hessian matrix given by 


- (Cf D n G^C 2 ) - (Cl D n G^C 2 ) T (76) 

(C;C,) 0) = (-D lp G^C 2 ) T (C„ - D lp GC 2 ) 

+ (C p -D lp GC 2 ) T (-D lp G™C 2 ) (77) 


Computation of H a 

Similarly, the derivative of the homotopy function with respect 
to the homotopy parameter a is the N x 1 vector 



and using Eq. (53), 

W 2 = vec (^"{ 2 [(^D 12 GC 2 - DlA - B t 2 Q x )LCl 

+ ( kDl p D lp GC 2 - XD\ p C p - B\ L p )XCl ] (66) 

= vec \2[(Dl 2 b x2 G^C 2 - flj Q^)LC 2 
+ ( D T n b i2 GC 2 - Df 2 C, - Bl Q x )L u) Cl 


7- an 

- — - (78) 

3a 

md 

3 H f3(3£/3G)1 

da L 3a j 

= vec(j-{2[(bl 2 b l2 GC 2 - D[ 2 C, - B[ Q x )LC t 2 

+ (A .Dl p b lp GC 2 - kDl p C p - Bl L p )XC t 2 }]^ 

= vec {2[(f) 13 DijCC 2 + b\ 2 D l2 GC 2 — D , 2 Ci 

- Bl Q x )LCl + (bl 2 b l2 GC 2 - b T n c , - Bl Q x )tcl 
+ (kbl p D lp GC 2 - kDl p C p - Bl L p )XCl] } (79) 


+ (kbl p D lp G (j) C 2 - Bl Lf)XCl 
+ ( kbl p b lp GC 2 - kb[ p c p - Bl L p )X^Cl]j (67) 


(*) = 


£(*) 

da 


( 80 ) 
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Implicit in these equations is the assumption that all system matrices 
are fixed with the exception of D i2 and D\ p and the parameters X 
and y. In general, the homotopy can be performed with any/all 
system matrices deformed, but the methodology for mixed H 2 /H 00 
employed is to vary only one parameter at a time. 

The derivative terms in the third equality in Eq. (79) depend on 
the deformation undertaken in the specified problem, i.e., the initial 
and final problem. Suppose that the deformation of the matrix b l2 
is prescribed to be 

&12( a ) = D]2,o( a ) + 0t[Di2,/(0() — £>12.oW] ( 81 ) 

where the 0 and / subscripts indicate the initial and final system 
matrices, respectively. It follows that 

^12 = &12 ./ ” D 12.0 (82) 

The derivative of other plant matrices and the parameters X and y are 
determined accordingly. To obtain expressions for the derivatives of 
L, Qoc, L p , and X with respect to or, Eqs. (49-52) are differentiated, 
resulting in the following: 

0 = (a + y~ 2 BB T Q x )L + L(A + y- 2 BB T Q x ) T 
+ (rt + Lr r ) (83) 

o = (A + y- 2 BB T Q x ) T Q x + Q m (A + y~ 2 BB T Q x ) 

-( 2y~ 3 yQ x BB T Q x ) (84) 

0 = A t L p + L p A + (A T L P + L p A + * C T p C p 

4- XC T p t. p + ACjCp) (85) 

0 = AX + XA t + (AX + X A T + B p B T p + B p B T p ) (86) 

where 

T =X - 2y~ 3 y BB t Q x + y ~ 2 (BB t Q x +BB T Q x +BB r Q x ) 



(87) 

and, from Eqs. (45-47), 


A = A — B 2 GC 2 — B 2 GC 2 

(88) 

B p = *p 

(89) 


(90) 

C p = t p - 6 lp GC 2 - b Xp Gtz (91) 

t = t x - 6 l2 GC 2 - D n GC 2 (92) 

The use of the canonical compensator formulation in the aug- 
mented system dynamics simplifies not only the necessary condi- 
tions but also the derivative expressions of the system matrices. 
From Eqs. (37-40), the augmented matrix derivatives reduce to 

, r A f -Ao O' 

|_ — C 2 . 0 ) 0_ 

(93) 

j r B p f - Bp* 
p l-N°(D 2 p , f -D 2 p .o)_ 

(94) 

s r &i.f - #i.o 

1 L- Af0 ( I) 2i./-r> 21 . 0 )_ 

(95) 

g r ^ 2 ./ — # 2 ,o 0 

2 “ L -N\D n .f - D 22 . 0 ) 0 

'] < 96 > 

G p = [C Pr f — C Pi o 0] 

(97) 


C\ = [C u — C 1>0 0] 

(98) 

<?2 = [0 0] 

(99) 

b\ p = W\pj “ A/ 7.0 o] 

(100) 

b\ 2 = [D\ 2 f — D\ 2 q 0] 

(101) 


The presence of the zero subblocks significantly enhances the 
computational efficiency of this approach. When implementing the 
procedure described at the beginning of this section, the preceding 
equations may be further specialized. In the initial H 2 homotopy 
procedure the initial and final plant matrices are the same and the 
homotopy is performed only on the measurement and process noise 
intensities, D\ 2 and D 2 \- Hence A, & 2t C\> and C 2 are identically 
zero. For the mixed H 2 /H » homotopy, the and H 2 homotopies 

are performed distinctly, which simplifies the computations signif- 
icantly because the plant matrices remain fixed and only y or X are 
varied at one time. 

Design Example 

To demonstrate the homotopy algorithm applied to optimal con- 
troller synthesis, the four disk example originally described in 
Ref. 19 and more recently by numerous others 10 will be used. The 
four disk model used in the example problem was derived from a 
laboratory experiment and represents an apparatus developed for 
testing of pointing control systems for flexible space structures with 
noncollocated sensors and actuators. Four disks are rigidly attached 
to a flexible axial shaft with control torque applied to selected disks 
and the angular displacement of selected disks measured. The plant 
parameters are taken from Ref. 4. 

H 2 Case 

To demonstrate the homotopy algorithm applied to H 2 controller 
synthesis, a direct comparison between the homotopy algorithm of 
Ref. 10 called HAO and the FOCUS algorithm (called H2HOM 
for the H 2 case) will be presented first. The main distinction be- 
tween the two homotopy algorithms for the H 2 case is the compen- 
sator architecture. HAO employs a general architecture that may be 
restricted to various parametrizations including the controllability 
canonical form, which is similar to the controller canonical form 
used in H2HOM for this SISO example problem. The HAO code 
has been highly optimized for efficient computation with the result 
that superfluous computations are not evaluated. The homotopy al- 
gorithm of this paper is patterned after the general approach of HAO 
and utilizes some of the more efficient computational aspects of the 
HAO code. 

The control design philosophy for this example is to scale the 
nominal control weight and the nominal sensor noise intensity by the 
parameter q . As q is reduced, the control authority is increased. For 
comparison with the published results in Ref. 10, a full- (eighth-) 
order compensator was synthesized. Although the results can be 
directly obtained from the LQG Riccati equations, the full-order 
compensator was chosen to tax the H2HOM algorithm, which must 
optimize over a greater number of parameters with increasing com- 
pensator order. 

Table 1 shows a comparison of the results from the H2HOM and 
HAO algorithms for the full-order compensator along with results 
from the H2HOM algorithm for sixth- and second-order compen- 
sators. All pertinent parameters as well as logic for step-size scal- 
ing and the computation of the prediction and correction errors are 
identical in both algorithms, which are implemented in MATLAB 
on a 486 66-MHz computer. Whereas the HAO code required a 
minimum step size of 1.907e-7, the H2HOM solution was much 
better conditioned and required a minimum step size of 0.025. As 
a consequence of the smaller step sizes with HAO, 2504 Hessian 
computations were required, as opposed to only 63 Hessian com- 
putations with H2HOM. The HAO code has been tuned extensively 
for efficient computation, as is reflected in the small number of 
floating point operations per second (flops) required. In spite of the 
significantly smaller number of flops with HAO, the H2HOM code 
required significantly less clock time for convergence to the same 
final compensator. (The results generated by the authors using the 
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Ifeble 1 Comparison of H2HOM and HAO algorithms 


Algorithm 

HAO 

H2HOM 

H2HOM 

H2HOM 

Compensator order 

8 

8 

6 

2 

Number of Hessian 

2504 

63 

60 

30 

computations 

Min. step size 

1.907e-7 

0.025 

0.05 

0.1 

Max. step size 

0.1 

0.1 

0.1 

0.2 

Max. number of correction 

9 

7 

8 

5 

iterations 

Mflops 

287 

936 

455 

37 

Time, s 

5104 

883 

488 

73 


Table 2 Comparison of FOCUS and Ref. 4 results 

Y 

IIT ztt p lloc.BH 

lloo.f 

WTipwp 11 2, HH 


2.0 

1.18 



0.382 

— 

1.5 

1.06 

— 

0.389 

— 

1.0 

0.855 

0.947 

0.410 

0.398 

0.9 

0.797 

0.889 

0.420 

0.403 

0.8 

0.732 

0.776 

0.432 

0.421 

0.7 

0.661 

0.668 

0.451 

0.438 

0.6 



0.589 

— 

0.460 

0.52 

0.511 

0.520 

0.512 

0.508 

0.5 

— 

05 

— 

0.518 

0.4 



0.398 

— 

0.572 

0.3 

— 

0.300 

— 

0.799 


HAO code differ slightly from the published results, 10 although the 
parameters in the HAO algorithm are the same. It is likely that the 
published results were generated with an earlier version of the HAO 
code. The qualitative trends remain the same.) 

In Ref. 10, the controllability canonical form is assessed as poorly 
conditioned because of the small minimum step size. However, 
Table 1 indicates that the static gain formulation in H2HOM yields 
a substantial improvement in conditioning along the homotopy path 
over the HAO implementation of the canonical compensator. Al- 
though the static gain formulation is better conditioned for this ex- 
ample problem, this may not be the case in general. An even more 
significant benefit of this formulation is the straightforward exten- 
sion to the Hoc and mixed H 2 /Hoq problems. 

Mixed H 2 /H 00 Case 
One-Input/Two-Output Case 

The seminal paper 4 dealing with mixed H 2 /H 0 0 design addressed 
the case where u; = w p in Eq. (37) with results from the four disk 
problem given for the full-order case. In this section, the one- 
input/two-output case will be repeated with the FOCUS algorithm as 
well as a two-input/two-output case with full- and fixed-order com- 
pensators to demonstrate the capabilities of the homotopy algorithm. 

Table 2 presents a comparison of the results from the FOCUS 
algorithm and the results published in Ref. 4. In Table 2, the BH 
subscript indicates results from Bernstein and Haddad 4 and the F 
subscript indicates results from FOCUS. Gaps in the columns of 
Table 2 denoted BH correspond to values of y for which results 
were not published. The absence of data generated by FOCUS for 
y = 2 and 1.5 is due to the fact that y is larger than the maximum 
00 norm, as described in the next section. A key distinction between 
the formulations of the mixed H 2 j optimization problem used 
in this paper and that of Ref. 4 is that whereas their formulation 
minimizes an overbound on the H 2 norm, the formulation utilized 
in FOCUS minimizes the actual H 2 norm. Consequently, as shown 
in Fig. 2, the FOCUS algorithm results in smaller H 2 norms for a 
given y with the gap in the 00 norm overbound smaller than the 
Ref. 4 results in the meaningful region for y, as described in the 
next section. Whereas y = 0.52 was the smallest y value reported 
in Ref. 4, y values of 0.5, 0.4, and 0.3 are utilized with FOCUS, as 
indicated in Table 2. This example also demonstrates the synthesis 
of fixed-order mixed controllers with singular constraints. In 
practice, compensators representing the extreme values of y or X 
typically will not be used because either performance or robustness 
would be severely diminished. A compromise value in the elbow 



of Fig. 2 typically would be chosen. Although the FOCUS results 
in Table 2 are not significantly better than the results of Ref. 4, the 
example demonstrates that the FOCUS code performs satisfactorily 
and lessens the gap between the overbound and the 00 norm. 

Discussion 

This example serves to illustrate some interesting features of the 
mixed H 2 /Hoo formulation implemented in FOCUS, as well as dis- 
tinctions from other formulations. The formulation implemented in 
FOCUS is a method for generating suboptimal controllers of 
fixed order that are subject to an H 2 constraint The cost functional 
for the mixed H 2 fHoo problem can be written as 

Jmix = 7oo + kJ 2 (102) 

where 

Joo = Vc{Q x BB t \ (103) 

y 2 = tr [XClC p } (104) 

subject to the corresponding Riccati and Lyapunov equations. The 
resulting Lagrangian, also given by Eq. (48), is 

C = ti{Q 00 BB T +(A T Q 00 + Q 00 A +C T C + y^Q^BB 7 Q X )L 

+ \XC T p C p +(AX + XA T + B p B T p )L p ) (105) 

Minimization of Joo results in an H& controller with an upper 
bound on the 00 norm given by the 00 norm of the H 2 controller (as 
y — ► 00 , the H 2 compensator for T zw is recovered by minimizing 
Joo). By successively lowering gamma, the minimum H « norm 
controller for 7^ is obtained. Minimization of J 2 results in the 
optimal H 2 compensator for T zpw . Thus, when nonzero X is used in 
/ mix , the Hoc cost functional imposes an additional constraint on the 
H 2 norm of T zpw , and for large y, the necessary conditions for the 
mixed problem yield the simultaneous solution of two H 2 problems. 
By increasing the H 2 weight, X, for a fixed y, ||T V J| 2 is reduced 
while \\T zw \\oc approaches the gamma overbound. At that point, the 
minimum H 2 controller for r vu) such that Halloo < y is obtained. 

For this example problem, the optimal H 2 controller for T zpu) 
results in \\T zw \\oc = 1.392 and \\T ZW \\ 2 = 0.3786. The y values 
larger than 1.392 are not meaningful for our formulation because 
the optimal U 2 compensator for T zpw satisfies the y constraint on 
HTiJU In that case, the //«, norm constraint is inactive in the 
mixed norm optimization. However, since the Ref. 4 formulation 
seeks to minimize the H 2 norm overbound, it is possible for their 
formulation to generate meaningful solutions for y > 1 .392. When 
using FOCUS, the first step should be to establish an upper bound on 
y by computing ||r zu , ||oo resulting from the optimal H 2 compensator 
for Ttpv. 
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Table 3 Comparison of two-input/two-output results 


y 

II T ZU1 Hoc, 8 

II T zpwp H 2.8 

^Tzw Hoc . 4 

\\Tzpwp II 2 . 

2 





2 

0.421 

1.5 

1.304 

0.446 

1.483 

0.437 

1.0 

1.000 

0500 

1.000 

0.496 

0.9 

0.900 

0.529 

0.900 

0532 

0.8 

0.789 

0560 

0.800 

0.587 

0.7 

0.694 

0.652 

0.700 

0.640 

0.6 

0597 

0.799 

0.600 

0.952 

0.5 

0500 

3.842 

0500 

1.180 

0.4 

0.400 

8.281 

— 

— 

0.3 

0.300 

24.84 

— 

— 



Two- Input/Two- Output Reduced-Order Case 

To demonstrate the capability of the homotopy algorithm to syn- 
thesize fixed-order compensators for the two-input/two-output case, 
the four disk example was modified by creating a fictitious distur- 
bance input for the H^ minimization subproblem. The new coef- 
ficient matrices are given by B\ =2 x B p and D 2 \ — 2 x D 2p . 
Additionally, a small penalty was placed on the control in the H » 
subproblem (D n = [0 0.001] r ). Table 3 shows the results of the 
comparison for the full- and fixed-order two-input/two-output case. 
Figure 3 plots the intermediate values of Table 3 and illustrates the 
cost trade (extreme values are removed to preserve scale). In Table 3, 
the subscript 8 indicates the full-order design and the subscript 4 
indicates the fixed-order design. In general, full-order compensator 
synthesis is more taxing on the numerical algorithm due to the larger 
number of parameters to be optimized, which is evident in both 
Table 3 and Fig. 3. The fixed-order design is able to achieve smaller 
H 2 norms at the smaller y values. Also, the oo norms for the larger 
values of y are closer to the overbound with the fixed-order designs. 
The reduced number of parameters to be optimized results in more 
accurate solution of the minimization problem. The closeness of the 
Hi/ Hoo cost curves in Fig. 3 indicates that the fourth-order compen- 
sator closely approximates the performance/robustness tradeoff of 
an eighth-order compensator except at the low end of the attainable 
Hoo norm values. The infinity norm of T zw resulting from closing the 
loop with the full-order optimal H 2 controller for T zpwp is 1.845, so 
that the full-order case with y = 2 was not a meaningful solution. 
However, a meaningful solution for y — 2 was achieved with the 
fixed-order compensator. 

Conclusions 

A homotopy algorithm is developed to synthesize fixed-order 
H 2 , Hoc, and mixed H 2 /H< Xi compensators employing a controller 
canonical form, and a representative flexible structure is used to 
demonstrate the numerical results. Fixed-order mixed Hi/H^ com- 
pensators are synthesized for the case of one disturbance input set 


and two outputs sets, as well as the case of two disturbance input 
sets, two output sets. These numerical results indicate that for the 
example presented, the static output feedback formulation using the 
controller canonical form is a more efficient means of synthesiz- 
ing dynamic compensators than employing a general compensator 
architecture. The synthesized reduced-order compensators perform 
well when compared to full-order controllers. The fixed-order mixed 
H 2 j Hoo formulation is shown to offer improved performance over 
standard H compensators by minimizing the H 2 norm while re- 
moving (or reducing) the gap between the actual H ^ norm and the 
gamma overbound. 
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